X=load("out_gamma_099");s1=X(:,1)==1;s2=X(:,1)==2;s0=X(:,1)==0;  s3=X(:,1)==3;s4=X(:,1)==4;s6=X(:,1)==6; s8=X(:,1)==8; s10=X(:,1)==10; s12=X(:,1)==12; s14=X(:,1)==14;s11=X(:,1)==11; s13=X(:,1)==13; plot(X(s0, 5), ';serial;', X(s1,5), ';thompson;', X(s2, 5), ';mean;', X(s8, 5), ";HP;", X(s12, 5), ";mHP;", X(s10,5), ";BB;", 125*ones(16,1), ";baseline;"); legend("location", "northwest")